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Abstract 

In this paper, the study of bivariate generahsed beta type I and II distributions is 
extended to the complex matrix variate case, for which the corresponding density 
functions are found. In addition, for complex bimatrix variate beta type I dis- 
tributions, several basic properties, including the joint eigenvalue density and the 
maximum eigenvalue distribution, are studied. 

1 Introduction 

Complex matrix variate distributions play an important role in various fields of research. 
Applications of complex random matrices can be found in multiple time series analysis, 
nuclear physics, complex multivariate linear model, shape theory, the evaluation of the 
capacity of multiple-input multiple-output (MIMO ) wir eless communication syste ms, see 
Mehtal (|l99ll ). iKhatril (|l965h . iMicheas et ai] (|2006l ) and iRatnaraiah et all (|2005al ) among 
many others. 

Several studies on the distribution of complex random matrices have been made, see 
Mathai and ProvostI (2005). The complex matrix variate Gaus sian distribution was intro- 
duced bv lWooding (|l9563 . and further develop ed by Turin! ( 196d) and lGoodman ( 1963[ ). The 
complex Wis hart distribut ion was studied by GoodmaiT ( 19631 ) and James among 
many others. Jamed ( 19641 ) and Khatri ( 1965 ) derived both complex central and noncentral 
matrix variate beta distributions. 

When the goal is to generalise the distribution of a random variable to the multivariate 
case, two options are normally addressed, by which it is extended to either the vectorial 
or the matrix case, e.g. normal, t or bessel distributions, among many others. However, 
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some of these generalisations have traditionally been made directly to the matrix case, when 
such a matrix is symmetric, as in the case of the chi-square and beta distributions, for 
which the corresponding multivariate distributions are the Wishart and matrix variate beta 
distributions, respectively. Nevertheless, these latter generalisations are inappropriate in 
some cases, because, in some a pplications, the researcher is interested in a vector variate, 
not in a symmetric matrix, see iLibbv and Novick ] (|l982h . In other words, the researcher 
is interested in a vector, say, X — {x\, . . . ,Xm)' ■, such tha t Xj has a margina l beta, type 
I or II distribution for all « = 1, . . . , to. In this respect, iLibbv and Novick I (jl982l ) and 



Chen and Novick I (|1984| ) proposed a multivariate vector- version of the beta type I a nd II 
distributions. Let us consider the following bivariate version, see lOlkin and Liul (j2003l ). 

Let Xq, Xi, X2 be distributed as independent gamma random variates with parameters 
a = ao,ai,a2, respectively (see Definition 12.11 in Section 2); and define 



Xi 



Xi + Xq ' 



U2 



X. 



X2 + Xq 



(1) 



Clearly, Ui and U2 each have a beta type I distribution, Ui ^ S/i(ai,ao) and U2 ^ 
BIi{a2,ao), over < ui,U2 < 1 (see Subsection 2.1). However, they are correlated such 
that (C/i,J72)' has a bivariate generalised beta type I distribution over < ui,U2 < 1. The 
kernel of the joint density function of Ui and U2 is 



oc 



^(1 



Ul 



102+00-1 



(1 - UiU2)°i+°2+ao 



< Ul, U2 < 1. 



A similar result is obtained in the case of beta type II. Here it defines 



F2 = ^. 

Xo 



Once again it is evident that Fi and F2 each have a beta type II distribution, Fi ~ 
BIIi{ai,ao) and F2 ^ S//i(a2,ao), over /i,/2 > 0. As in the beta type I case, they 
are correlated such that {Fi,F2y has a bivariate generalised beta type II distribution over 
/i, /2 > 0. The kernel of the joint density function of Fi and F2 is 



/ai — 1 ra 
1 J2 



(1 + /l + /2)"i+''^+"" ' 



./l,/2>0. 



Some appli cations to utility modellin g and Bayesian analysis are presented in lLibbv and Novick 
(|l982l ) and lChen and Novick I (|l984l ). respectively. 

These ideas can be extended to the matrix variate case. Thus, let us assume a parti- 
tioned matrix U = (Ui:U2)' G C^mxm^ then under the complex matrix variate versions of 
the transformations (H]), we are interested in finding the joint density of Ui and U2, from 
where it is easy to see that the marginal densities of Ui and U2 are complex matrix variate 
beta type I distributions. In the central and noncentral real cases, the matrix variate joint 
densit ies of Ui and U2 and of Fi and F2 , togeth er with some of their properties, are stud- 
ied in Dfaz-Garcia and Gutierrez- JaimezI (|2009allbl ) . These distributions are termed central 
complex bimatrix variate generalised beta type I and II distributions, respectively. 

In the present paper, the bivariate generalised beta type I and II distributions are ex- 
tended the the complex matrix variate case, see Sections [3] and 01 In Section [51 certain basic 
properties, the joint eigenvalue density and the density of their maximum eigenvalues are 
studied for the complex bimatrix variate generalised beta type I distribution. 
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2 Preliminary results 



The join eigenvalues density can be calculated using complex hypergeometric functions and 
complex invariant polynomials with matrix arguments. In this section, these and other issues 
are addressed for the case of the complex multivariate distribution. Let us first establish 
some notation. 



2.1 Notation and matrix variate distributions 

Let A — {ars) — {arsi + iars2) — -^i + be an m x rt of complex numbers, A e (rmxn^ 
where Ai — {ars^) and A2 = (a^sa) are real matrices mx n, Ai, A2 G 3?™^", and i = 
Then, A' denotes the transpose of A, A denotes the conjugate of A, and A^ denotes the 
conjugate transpose of A. For n — m, let tr(A) = an + • • ■ + amm, etr(tr(A)) = exp(tr(A)), 
and then |A| denotes the determinant of A. A = A^ > is a Hermitian positive definite 
matrix, and A^/^ denotes the unique Hermitian positive definite square root matrix of 
A = A^ > 0. dA denotes the differential matrix of A and {dA) = {dAi){dA2) denotes 
the volume element (Lebesgue measure) commonly associated with A. For example if A is 
a Hermitian matrix, then Ai is a symmetric matrix (A = A') and A2 is a skew-symmetric 
matrix (A = — A'). In this case 

(dA) = f\ darsi /\ dars2 ■ 

r<s r<s 

The space of all matrices Gi e (m < n) with orthonormal columns is termed the 

Stiefel manifold, denoted by £V„i.„. Thus 

e:v™,n - {Gi e IGfGi = i„,}. 

From ljames (1963) 

Vo1(g:V™,„)= / (GfdGi) = 



/GiG£V,„,„ CTmN' 

where CFm [a] denotes the complex multivariate gamma function and is defined as 

£r„[a] = / etr(-V)|V|"-"(dV) = 7r™("-i)/2 TT r[a - j + 1], 

Jv=v«>o 

where Re(a) > m — 1. 

If m = n, it is have a special case of the Stiefel manifold, termed unitary manifold or 
unitary group and denoted as U{m) = €Vm,m.- 

Definition 2.1 (The complex multivariate beta function). The complex multivariate beta 
function, denoted as £/3„Ja, 6], is defined by 

e:/3„j6,a] = / |sr-™|i„,-s|^-™(ds) 



g:F„,Mg:F„j6] 

€^„^[a + b] 
where Re(a) > m — 1 and Re(6) > m — 1. 



/ |Rr'™|I™ + Rr('^+'')((iR), R = (I - sy 
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We now give definitions of the compl ex m atrix var i ate ga mma, beta type I and II distri- 
butions, see ljamed (|l964l ). iKhatoil (|l965l) andlMathaJ (|l997l) . 

Deflnition 2.2 (Complex matrix variate gamma distribution). It is said that A e (rrnxm^ 
random Hermitian positive definite, has a complex matrix variate gamma distribution with 
parameters a and a Hermitian positive definite matrix G (rmxm^ j£ j^-g (jgj^gjty function is 

^j-ij^|Ar™etr(-e-A)(dA), A = A- > 0, (2) 

where Re(a) > m — 1. Such a distribution is denoted as A €Gmia, ©)■ 

Lemma 2.1 (Complex matrix variate beta type I and II distribution). If A and B have 
a complex matrix variate gamma distribution, i.e. A ~ £^,„(a, !,„) and B ^ €Gmib,lm) 
independently. 

1. Then the complex matrix variate beta type I distribution is defined as 

( (A + B)-i/2a((A + B)-i/2)', Definition I or, , . 

\ Ai/2(A + B)-i(Ai/2)', Defimtion2. ^ ' 

Thus under definitions 1 and 2 its density function is denoted as 

€BI,niU;a,b), 

and given by 

^ |Ur-™|I„-U|''-™(dU), 0<U = U«<I,„, (4) 



(5) 



e:/3^[a,5]' 

this being denoted as \J ^ <LBIm{a, b) with Re{a) > m — 1 and Re{b) > m — 1. 

2. Then the complex matrix variate beta type II distribution is defined as 

B-i/2A(B-i/2)', Definition 1, 
Ai/2B-i(Ai/2)', Definition 2, 

Thus under definitions 1 and 2 its density function is denoted as 

€BII,n{V;a,b), 

and given by 

c^R ] ,J Fr-'"|I,„+ F|-(°+'').(dF), F = F^>0. (6) 
this being denoted as F ^ €BIIm{a, b) with Re{a) > to — 1 and Re{b) > m — 1. 

2.2 Complex hypergeometric functions and invariant polynomials 

The following defin ition s of hypergeometric functio ns with a matrix argument are based on 
Constantinj (|l963h and iKoev and Edelmanl (|2006f) . 



Definition 2.3. The hypergeometric functions of a matrix argument are given by 

pFl^ ){ai,...,ap;bi,...,b,:X) = ^}_^-—^ — ^ , (7) 

t=0 T {Ol)r ■■■[Oq)r 
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where denotes the summation over all the partitions r = (ii, . . . , t„i), ti > ■ ■ ■ > tm > 0, 



of t = ii + • • • + trm Ct"''(X) is the Jack polynomial of X corresponding to t and the 
generalised hypergeometric coefficient (a)?"'* is given by 



where (a)t = a(a + l)(a + 2) • • • (a + <— 1), (a)o = f . Here X S is a complex symmetric 

matrix and the parameters a^, 6^ are arbitrary complex numbers. 

Other charac teristics of the parameters aj and b j an d the convergence o f (l7l) a. ppear in 
MuirheadI (|l982l . p. 258). lGross and Richardl (|l987l ) and lRatnaraiah a/.l (|2005bl ). 



Remark 2.1. In Definition 12.31 when a = 1 and 2 the complex and real cases are ob- 
tained, respectively. In this paper it is considered only the complex case. Then, adopting 
the notation used by I Jamei (1964), denoting the complex hypergeometric function, gener- 
alised hypergeometric coefficient and zonal polynomials as pFq"^ = pFq^\ [a]r = (a)r^', and 
Ct(-) = Ct^\-), respectively 

A special case of ([7]) is 

iFo(a;X) = ^^[a].^ip^, (||X|| < 1) 

t = T 

= |I™-xr^ 

where || X|| denotes the maximum of the absolute values of the eigenvalues of X. 
From Ratnarajah et al. ( 2005lJ l. is known that. 



ii^i(a;c;X) = 



1 



X ( o-Fo(XY)|Y|''-"|I- Y|"-"-™(dY), 

"'0<Y=Y^<I„ 



and 



2Fi(a, ai;c; X) 



^/3m[a,c- a] 

X / i^i)(ai;XY)|Y|''-™|I- Y|'=-''-'"(dY). 

JO<Y=Y"<I™ 



(8) 



Thus we have 



Lemma 2.2. Let X = X^ e g:™x", with Rc(X) < I, Rc(a) > m-l, Re(c) > m-l and 
Re(c — a) > m — 1. Then 



-iFq+i (a, ai, • ■ • , Op; c,bi, . . . ,bq; X) 



•^/^mfa^c- a] 

pFqiai • • • ap; 6i • • • 5,; XY)|Yr-'"|I - Y|'=-''-'"(dY). 

/0<Y=Y«<I,„ 

Proof. First, an expansion is applied in terms of complex zonal polynomials 
pF,(ai,...,ap,6i,...,5„XY)-^2. • • • [fe,]. t\ ' 
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Then, after integrating term by term, see lRatnaraiah et al\ (|2005br ). we have that 
pF,(ai • • • ap; 61 • • • 6,; XY)|Yr-™|I - Y|^-'^-™(dY) 



0<Y=Y«<I„ 



ai 



|Y|a-m|I_ Y| 



rt-/^ r 1 \ " \ " Mr [Ql]r ■ ■ ■ [Qp]r ^ 



c—a—m 



Cr(XY)(dY) 



t=0 r 



= £/3,„[a, c - a] p+ii^g+i(a, ai • • • Op; c, 5i • • • 6^; X), 

and the required result follows. □ 
The use of complex zonal polynomials and the hypergeometric function with a matrix 
ar gument has only rece ntly been exten ded; to a lar ge extent this is derived from the work 
of Koev and Edelman 

] 12006), who in iKoevI (|2004l) provided a program in MatLab with 
a very efficient algorithm for calculating Jack polynomials (in particular complex zonal 
polynomials) and the complex hypergeometric function with a matrix argument. 

This section concludes by establishing the following two properties of a class of homoge- 
neous polynomials C^'^(R, S) of degrees k and t in the eigenvalues of the Hermitian matrices 
R, S e e:'"^™, respectively, see lDavisI (|l980f ) and lRatnaraiah et a l} (l2005 a||. Th ese proper- 
ties generalise the in c ompl ete beta function equation (61) of IConst antina (|l963l ). The first 
is proposed by David ( 19791 . eq. (3.3)) and the second is obtained using the complex versions 
of Chikuse" (Il980l eg. (3.3 3)) and the review version of IChikuse ( 19801 . eq. (3.11)) given in 



[Chikusc an d David ((l986l eq. (2.7)) 
Lemma 2.3. Let R, S, fJ, H e (rmxm ffg^rmitian matrices. Then 

,A=A" 



/ 

Jo 



IRI' 



II -R 



b—7n 



a(f^R)(dR) 



A;— K\(I>^K..T 



k\[a + m](j, 



(9) 



L 



|R|°-"|I - R|''-"C^'^(H, nR)(dR) 

H^, [-b + m].<;'-^^^c;;'=--(A, H, nA) 



(10) 



sl[a + m]0j 

where, from Chikuse and Davik 1 198t . Lemma 2.2(i) and (ii), respectively) it is known that 
5^'^ (I, I) 



3 Bimatrix variate generalised beta type I distribution 

Let A, B and C be independent complex random matrices, such that A ~ €Gm{(i,im), 
B ~ (ZQm{b, Im) and C ~ £fj,„(c, I,„) with Re(a) > to— 1, Re(6) > to— 1 and Re(c) > to— 1 
and let us define 



Ui = (A + C)-i/2a(a + C)-i/2 and U2 = (B + C)-i/2b(b c)-i 



/2 



(11) 
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Of course, Ui ^ CBImia, c) and U2 ^ €BI„i{b, c). However, they are correlated such that 

the distribution of U = (Ui:U2)' G (f'^ni'xm ^g^j^ ]-,g termed a complex bimatrix variate 
generalised beta type I distribution, denoted as 

V ^itBgBl2,nxm{a,b, c). 

Theorem 3.1. Assume that U ^ BgBl2mxm{o-, b, c). Then its density function is 

iuii-™|U2i^-"Mi,„ - Uiii^+^^™|i™ - U2r+^-™ 



e:/3:„[a,6,c]|I„-UiU2h+''+^ 
< Ui = < Ijn, < U2 = < Im, where the measure 

(dU) = (dUi) A (dU2), 

and 

r 1 e:r„,[a]e:r„,[6]G:r„,[c] 

€f3,^[a,b,c\ = 



-(dU), (12) 



eTm[a + b + c] 
and Re(a) > m — 1, Re(6) > m — 1 and Re(c) > m — 1. 

Proof. The joint density of A, B and C is 

|^|Q-m|g|fc-m|Q|c-m 



£r,„[a]£r,„[6]G:r 

7) 



■ etr(-(A + B + C))(dA)(dB)(dC). 



By effecting the change of variable (llip , and taking into account iMathail (|l997l Theorems 
3.5 and 3.8, pp. 183 and 190, respectively) we have 

{dA){dB){dC) = |C|2™|I,„ - Uir^^ll,,, - U2|-2-(dUi)(rfU2)(dC). 

The joint density of Ui, U2 and C is 

|Ul|°~'"|U2|''~"' |p,|a+b+c-m 

X etr [-(I™ - U2)-i(I™ - UiU2)(I™ - V.y^C] {dC){dV,){dU2). 
Integrating with respect to C using 



r Icr+^'+^-'-etr [-(I„, - U2)-i(I™ - UiU2)(I™ - lJi)-^C] (dC) 

= €T[a + b + c] 



'C=C«>0 

ii,„-uir+''+=ii™-u2r+^+^ 



|i™-UiU2h+^+- 

(from ([2])) gives the st ated marginal density function for ( Ui :XJn )'. □ 
As in the real case ( Diaz-Garcia and Gutierrez- JaimezI l2009a ) . the joint density (fT2|) can 



be represented as a mixture. Let us first note that 

|I,„-UiU2r('^+''+^' = iFo(a + & + c;UiU2) 

= !Lz^[a + ^ + c]r 

t=0 r 
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By substituting in it is obtained that the joint density function of (Ui:U2)' is 

a+c-™ ^.(UlU^) 



Moreover 



x'|i,„-u2r+^ 



t! 



(13) 



^^[a + b + c]r€Tm[b + c]€rm[a + c] m u, \ ~ , . 



4 Bimatrix variate generalised beta type II distribution 

Let A, B and C S (rmxm independent, where A ~ £C/m(a,Im), B €^m(6, I„i) and 
C ^ €Qm{c, I„i) with Re(a) > to — 1, Re(6) > m — 1 and Re(c) > to — 1 and let us define 

Fi = C-1/2AC-1/2 and F2 = C-^^^BC-^/^. (14) 

Clearly, Fi ^ €BIIm{a,c) and F2 ^ €BIIm{b,c). But they are correlated and so the 

distribution of F = (Fi:F2)' S S)fj2mxm ^^^^ termed a complex bimatrix variate generalised 
beta type II distribution, which is denoted as F ~ €BQBI l2m-Km{o,, b, c). 

Theorem 4.1. Assume that F ~ BQBl2m,y.m{o,, b, c). Then its density function is 

-m, (15) 



iFi^-^lFal^-' 



€/3:ja,6,c]|I„ + Fi+F2h+^+^ 
Fi = > 0, F2 = ¥2 > 0, where the measure 

id¥) = (dFi) A (dF2). 
and Re(a) > m — 1, Re(&) > m — 1 and Re(c) > m — 1. 

Proof. As an alternative to proceeding a s in Theoreml3.1[ let us recall that if U ^ QiBIm (« , 6) , 

then ( Im— U)~^— Im ^ ^BIIm{a,h), see Srivastava and Khatri ( 1979t ) and Dfaz-Garcfa and Gutierrez- Jaime j 
(|20071 ). Then 

/Fi \ _/ (I,„-UO 

" V F2 ) ~ \ (I™ - U2) 

with the Jacobian given by fsee lMathail (Il997l . Theorem 3.8, p. 190)) 



(dUi)(dU2) 

Also, note that 

Im — (Ir?j + Fi) = 
Im — (Im + F2 ) = 



Fil 



-2m I 



-2m 



(I,„ + Fi)-i((I„,+Fi)-I 
(I„ + F2)-iF2, 



(dFi)(dF2) 

(Im -\ 



Fi)-iFi 



Then the joint density of (Fi:F2)' is 

|Fi|°"'"|F2|''-'"|Im+ri|-(°+^+°)|I^+F2|-(°+^+") 

G:/3:„[a, 6, c]|I„, - (I„, + Fi)-iFiF2(Im + ¥2)-^^+^+^ 
The desired results then follows, noting that 

|Im+Fi|-l|I,„+F2|-l , ^ , ^ ,-1 



(dF) 



□ 



|I„-(I,„ + Fi)-iFiF2(Im + F2)-i| ' 
Other properties of the distribution €BQBI l2my.m{0', b, c) can be found in a similar way. 



8 



5 Properties 

In this section, several basic properties and eigenvalue distributions are studied. 
5.1 Basic properties 

As direct consequences of Lemma 12.21 the following basic properties are derived: the mo- 
ments £'(|Ui|''|U2|^), the distributions of the product Z = U2^'^UiU2^^ and the inverse 



(Ur^iU ^M- Their proofs are similar to those given for the real case, see lDiaz-Garcia and Gutierrez- Jaimea 
()2009al) . 

Theorem 5.1. Assume that (Ui:U2) ^ <^BQBl2rnxrn{'^,b,c) then 

,, g^rn [a + r,b + c] €(3„-^ [b + s,a + c] 
-C'l Ui U2 j = — — 

X 3i^2(a + r, b + s , a + b + c; a + b + c + r, a + b + c + s; !,„), 
with Re(5 + r) > m — 1, and Re(a + c) > to — 1. 



Theorem 5.2. Consider that (Ui:U2) ~ €BGBl2mxm{0'ib, c). Then the density function 
of Z = Z" = Vl^^ViVl^^ e is 

€/3„Ja + c,6 + c]|Zr-"|I,„-Z|' 



-2^1 (a + c,a + c;a + b + 2c;Im~ Z)(dZ), 



Zr = ^^^21^ ' , ! 1^ —zF2{c, a + c,a + c;a + c + r,a + b + 2c- 1,„ , 

TOt/i < Re(Z) < Im, Re(a + 6) > m — 1 and Rc(& + c) > to — 1. 

Theorem 5.3. Let (Ui:U2)' ~ itBQBl2mxm{a,b,c). Then the density function of Y 



-(dV), 



ivir°'"iv2r'""'iim - vr'ii''+^-"'iim - v2"Y+^ 

G:/3:;,[a,6,c]|I™-(ViV2)-ih+''+^ 
< Vi = Vf^ < I„i, < V2 = \^2 < Inij where the measure 

{d¥) = (rfVi) A (rfVs). 
and Re(a) > m — 1, Re(6) > m — 1 and Re(c) > to — 1. 

5.2 Joint eigenvalue distribution 

Many statistics of the multivariate test hypothesis are functions of the eigenvalues or of the 
maximum eigenvalue. In these fin al two subsec tions, we find the eigenvalues distribution. 
The following result is needed, see JamesI ( 19641) . 
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Lemma 5.1. // /^(A) (dA) is the density function of a Hermitian matrix variate A € 
ffrnxm^ i/ien the distribution of the diagonal matrix A = diag(Ai, . . . , Am), Ai > • • • > Am > 
0, of the eigenvalues of A, where A — GAG''^, G S U{m), is the eigendecomposition of A, 
is 



fJA) = ——^ 1[{X,. - A,)2 / /,(GAG^)(dG), 

G-rmH Ju{m) 

where {dG) is the invariant measure on the unitary group lA{m) normalised, given as 

(rfG) = ^^"'[^^ (GdG^), such that j {dG) = 1. (16) 

2m^m JU(m) 

Theorem 5.4. Assume that (Ui:U2) ^ itBGBl2mxm{a,b,c) and let 



The spectral decomposition of Ui and U2, with E, G € U{m) and T)\ = diag(Ai, . . . , Am), 
1 > Ai > . . . , Am > and D5 = diag((5i, . . . , Sm), I > Si > . . . ,S„i > 0. Then the joint 
density function of Xi, . . . , Am , . . . , 6m is 



^2m(m— 1) 



r<s e<f k=0 k ' '~^K.(i-m) 

with Re(a) > m — I, Re(5) > to — 1 and Re(c) > 771 — 1. 

Proof. . The desired result follows b y making the t ransformation (|17p . from Lemma 15.11 
equation p6|) and observing that, see lChikusd (|2004l ) 

[ [ C,iEn^E"GT)sG")idE)idG) = ^'^(DA)g.(D,) ^ ^ 

5.3 Joint distribution of Xmax cind (Jmaa; 

In this subsection, we derive the distribution of the largest eigenvalues, Amaa; and 5max of a 
complex bimatrix variate beta type I matrix. With this aim, let us consider the following 
theorem. 

Theorem 5.5. Assume that (UiiUa) - €BgBl2mxm(a,b,c) and let Ai,A2 G 

be Hermitian positive definite matrices, < Ai, A2 < I. Then the probability P(Ui < 
Ai,U2 < A2) is given by 

€(3m[a,m]€f3m[b, m] V V Mfe M^i [a + b + c]^ 

tf3*Ja,b,c\ I 111 21 2^ 2^ ^, ^, ^, 

^ [-(a + C) + to]^[-(6 + c) + ?7l]^ <T,r.K;02p^<T,T,K. . A A A ^ 

" [a + TO],J5 + TO],, ^^^"-^^ (A„A2,A,A2) 

w/iere 

E E E 

anrf Re(a) > m — 1, Re(6) > tti — 1 anrf Re(c) > jn — 1. 
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Proof. From ^ the probability P(Ui < Ai,U2 < A2) is 



Ai 



o Jo 



lUil- 



lU, 



|b-r?j. I 



Ui||^+^-™|I, 



\a-\-c—m 



G:/3:„[a,6,c]|I™-UiU2h+^+^ 



-(dUi)(dU). 



(18) 



By rewritten (fT8|) as in (fT3|) . and integrating with respect to U2 using it is obtained 
that 

^(^m[b,m] .^^ ^ ^ [5]0ja + 5 + c]«[-(a + c)+m],0;f 



iuir-"ii,„-ui 



\b-\-c—my-iT,K 



O 



C;;''(A2,A2Ui)(dUi). 



The desired result follows by integrating with respect to Ui using ([T 
The following result is obtained from Theorem 15.51 



□ 



Corollary 5.1. Let (Ui:U2) ^ €.BQBl2mxm{a,b,c), Re(a) > m - 1, Re(6) > m - 1 and 
Re(c) > TO — 1. If Xmax Srnax the largest eigenvalues 0/ (Ui and U2, respectively, 
then their joint distribution function, P{\max < 2;, 5max < y) is given by 



€13 

rn 

€f3*J_a,b,c] 



E E 



/c! t\ s! 



-(a + c) + TO]^[-(b + c) + to]<^ ;,,^^^^^^K;02gK.T.<T;^CT,T,K 
[a + 6 + c]„[a + m]^, [6 + to]^, '""■^^ "^^ 



Proof. . The inequalities A,„Qa; < x and (5„ia2; < y are equivalent to Ui < xl and U2 < yl. 
Therefore, the result follows by lett ing Ai = xl and A2 — yl in Theorem 15. 5[ and taking 
into account that from fPavij ( 19791 . eqs. (2.1) and (2.7)) 



C^f^''(xl,yl,xyl) = x^+^y*+^-C^f''^(I). 



□ 



6 Conclusions 

Random matrices play a prominent role because of their deep mathematical structure. They 
have arisen in a number of fields (statistics, graph theory, stochastic linear algebra, physics, 
signal processing, etc.), often independently. In particular, the statisticians are interested in 
studying the Normal (Gaussian), Wishart, manova, and circular random matrices, which, 
from the point of view of r andom ma trix theory are termed He rmite, Laguerre, Jacobi, and 
Fourier ensembles, iMehtal (|l99ll) and lEdelman and Rad 

The results obtained in the present work can be considered as a generalisation of the 
Jacobi ensemble for the case in wh ich there are two correlated Jacobi ensembles, and with 
all their potential applications, see lEdelman and Sutton ( 2008 ). 



Another potential use appears in the context of complex shape theory, see lMicheas et al 



(l2006f). Specificallv Jji the approach known as afline shape or configuration densities, see 
Caro-Lopera et al\ (|2009l) . 
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